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ABSTRACT 

The  kinetics  of  free  radical  polymerization  were  treated  theoretically  under  the 
assumption  that  the  reaction  temperature  increases  at  constant  rate  (T  =  const).  From 
The  Arrhenius  relationships  for  the  kinetic  rate  constants,  approximate  corresponding 
functions  of  time  were  derived,  which  transformed  the  kinetic  equations  accordingly. 

A  differential  equation  describing  the  sum  of  the  concentrations  of  all  free  radicals 
appearing  in  such  a  system  was  solved  by  approximation,  and  from  the  resulting  series 
the  conditions  for  the  lower  as  well  as  the  upper  bound  of  the  stationary  state  assump¬ 
tions  were  obtained.  Defining  the  required  precision  by  a  large  number  N,  where  the 
relative  error  caused  by  the  stationary  state  assumptions  cannot  exceed  the  quantity  1/N, 
the  stationary  state  criterion  obtained  could  be  expressed  by  the  following  inequality: 

1 M  *t 

k<i 

where  f  is  the  efficiency  constant,  [  i],^  and  kt  are  the  functions  of  time  defining  the 

corresponding  concentration  of  initiator,  and  the  kinetic  rate  constants  for  initiator 
decomposition  and  termination. 

If  the  constant  rate  of  temperature  increase,  T,  approaches  zero,  the  reaction  becomes 
a  polymerization  carried  out  at  constant  temperature.  Under  such  conditions  the  derived 
criterion  turns  into  another  criterion  valid  for  constant  temperature.  This  new  expression 
is  identical  to  the  corresponding  formula  derived  previously  (References  1,  2,  and  3). 
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INTRODUCTION 

In  our  previous  work  (References  1,  2,  and  3)  a  stationary  state  criterion  for  free 
radical  polymerization  was  derived  from  the  rigorous  solution  of  the  kinetic  equations 
characterizing  the  system.  The  mathematical  analysis  was  carried  out  under  the 
assumption  that  a  constant  temperature  was  maintained  during  the  polymerization.  Since 
the  kinetic  rate  constants  are  functions  of  temperature,  according  to  the  Arrhenius 
relationship,  it  is  important  to  evaluate  the  appropriate  parameters,  by  carrying  out 
the  experiments  at  different  temperatures. 

During  the  last  decade  many  attempts  have  been  made  to  study  the  kinetics  of  free 
radical  polymerization  at  varying  temperatures,  and  to  derive  the  Arrhenius  parameters 
from  a  single  experiment.  It  was  found  that  the  most  favorable  technique  was  to  increase 
the  temperature  of  reaction  at  constant  rate  during  the  polymerization  (References  1 
through  5).  The  growing  importance  of  this  method  calls  for  the  extension  of  the  mathe¬ 
matical  analysis  made  previously  to  include  also  the  new  experimental  conditions.  Such 
a  treatment  was  carried  out  in  the  current  work  and  a  generalized  stationary  state 
criterion  was  derived.  The  new  formulas  obtained  are  valid  if  the  rate:  dT/dt  =  const, 
where  T  =  temperature,  and  t  =  time.  If  this  rate  approaches  zero,  the  new  expressions 
for  the  stationary  state  become  the  old  formulas  obtained  previously  (References  1,  2, 
and  3)  for  free  radical  polymerization  at  constant  temperature. 

KINETIC  SCHEME  FOR  FREE  RADICAL  INITIATED  POLYMERIZATION 

The  process  of  free  radical  polymerization  may  be  considered  to  involve  the  following 
steps  (Reference  6): 

Initiator  decomposition  and  initiation 


Propagation 


Termination 


R-  +  M 


•2R- 

k„ 


M-  +  M - — ►  M- 

n  n+ 1 


M-  +  M'  - M 

n  m  n+m 


M-  +  M-  - — — ►  M  +  M 

n  m  n  m 


R-  +.  R* 


t  r 


RR 


R*  +  M* 


'tn 


RM, 


where  kjjjka,  kp,  k^. . .  are  the  corresponding  kinetic  rate  constants. 


(1 ) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 
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Similar  to  the  method  utilized  previously  (References  1,  2,  and  3),  the  rate  constants 
of  termination  will  be  assumed  as  kxj.  =  and  k^n  =  2kj..  The  relatively  small  contri¬ 
bution  of  the  terms  ktr[R-]2andkta[R-][M-]  to  the  total  termination  makes  the  possible 
error  introduced  by  such  an  assumption  negligible.  The  system  defined  by  Equations  1 
through  7  leads  to  the  following  differential  equations: 


(8) 


(9) 


Adding  both  sides  of  Equations  8  and  9,  one  can  obtain: 

“57"  =  2kdM-2V2 

where 

r  =  [r  ]  *  M 

is  the  sum  of  concentrations  of  all  free  radicals. 

The  rate  of  the  initiator  decomposition  derived  from  Equation  1, 


leads  to  the  following  expression  for  [i  j : 

[i]  =  Io  e*p(-/  kddu) 

0 

where  Iq  is  the  initial  value  of  [i]  .  From  Equations  10  and  13  one  can  obtain: 


(10) 


(ii) 


(12) 


(13) 


t 

^- =  2kdexp  (-/  kddu)-2kfr2  (14) 

o 

Analogously  to  the  method  applied  in  the  previous  work  (References  1,  2,  and  3),  the 
stationary  state  will  occur  if  the  rate  dr/dt  in  Equation  14  is  much  smaller  than  either 
of  the  other  two  terms  appearing  in  this  equation.  With  such  an  assumption,  dr/dt  can 
be  omitted,  and  the  concentration  can  be  expressed  as  follows: 

1/2  t 

)  exp  (  y/  *<ddu  )  (15) 

o 
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Since  such  an  assumption  is  only  justified,  if  r  is  very  close  to  rg,  another  condition 
has  also  to  be  fulfilled  simultaneously: 

-r—  =  I  -  h 


1 1  -  -^1  =  Ihl  <  -i-  ( 

V  -  N 

where  N  is  a  large  number  defining  the  degree  of  precision,  and  l/N  is  the  maximum 
deviation  of  r/rg  from  unity.  The  criterion  for  the  stationary  state  can  be  derived 
according  to  the  conditions  defined  by  Equation  17,  if  the  expression  for  r  is  known 
explicitly.  To  obtain  such  an  expression  one  has  to  solve  differential  Equation  14. 

It  has  to  be  borne  in  mind,  that  the  temperature  of  the  system  under  consideration 
is  no  longer  constant,  and  must  be  defined  as  a  function  of  time  as  follows: 


T  =  T0  +  Tt 


where  T0  is  the  initial  temperature  and  T  =  (dT/dt)  =  const.  According  to  the  Arrhenius 
relationship,  any  kinetic  rate  kg  can  be  expressed  as: 


k.  =  A  exp 


(-fe) 


where  the  constant  parameters  Ag,  Eg,  and  R  are  respectively,  frequency  factor, 

activation  energy,  and  gas  constant.  T  is  expressed  in  °K.  Since  T  varies  with  time, 
kg  is  also  a  function  of  time.  This  function  defined  in  Equation  19  will  now  be  expanded 
into  the  following  series: 


In  general  (T  t/T0)«l,  and  Equation  20  reduces  to: 


Two  quantities,  constant  during  a  single  experiment,  will  now  be  introduced.  The 
initial  value  of  ks  at  T  =  T0: 


k  =  A  exp 
s,o  s 


(-ft) 


s  RT„ 
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Applying  these  two  experimental  constants,  the  Arrhenius  formula,  Equation  19,  can  be 
transformed  into  the  following  approximate  function  of  time: 

k  =  k  exp(f  t)  (24) 

S  S,0  s 

which  is  simple  and  convenient  for  further  calculations.  In  general,  the  term  (T  t/T0) 
appearing  in  Equation  20  is  a  small  number  over  the  entire  period  of  the  reaction,  and 
the  approximation,  Equation  21,  is  justified. 

Applying  the  formalism  given  by  Equations  22,  23,  and  24  one  can  express  the  kinetic 
rates  of  and  as  follows: 


where 


and 


where 


kd  =  kd,o  exP(fdf) 

•  (if) 

kt  =  \o  exP(ftt} 


(25) 


(26) 


(27) 


(28) 


(29) 


■.  •© 


(30) 


Equations  25  through  30  can  be  used  to  transform  Equation  14  into  a  new  equation: 


=  2k.  I  expf-^6-  +  f  t  -  )  -2k  eftf 

dt  d,o  o  ^  f  d  fd  /  fi° 


(3D 


Then  by  dividing  both  sides  of  Equation  31  by  k^,  and  introducing  a  new  variable,  x, 
one  gets 


(32) 
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and  the  following  differential  equation  is  obtained: 


dr 

dx 


ae 


-bx 


-P 


2 

r 


(33) 


where  p,  a,  and  b  are  constants  defined  as  follows: 


and 


(34) 

(35) 


(36) 


In  general  pel,  and  |  l-p]«l.  According  to  the  initial  conditions,  x0  -  kq.o/fq,  and 
r(xo)=0.  For  reasonable  values  of  T,  xQ  is  usually  a  small  number.  In  the  current  work 
it  will  often  be  considered  as  10"2. 

APPROXIMATE  SOLUTION  OF  DIFFERENTIAL  EQUATIONS 
OF  THE  SUM  OF  FREE  RADICALS 

Equation  33  was  solved  rigorously  in  the  previous  work  (References  1,  2,  and  3)  for 
p  =0.  For  p  1  O  no  rigorous  solution  was  achieved.  However,  methods  exist  which 
provide  an  approximate  solution  of  Equation  33,  and  enable  one  to  derive  a  criterion  for 
the  stationary  state;  one  of  these  methods  will  be  introduced  in  the  following: 

Under  stationary  state  assumptions,  Equation  33  reduces  to: 


“X  -p  2 
a©  -bx  r.  =  0 
s 


(37) 


Hence: 


(38) 


The  function  r  appearing  in  Equation  33  will  now  be  replaced  by  a  new  function,  g,  defined 
as 


(39) 


From  Equations  33,  38,  and  39  one  can  derive  a  differential  equation  for  g: 

do 

-j^-=  -2c<?  2  x2  g  +  i-  (-g--l)d-g2)  (40) 

l/2 

where  c  =  (ab)  . 
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The  relationship  given  by  Equation  39  can  also  be  expressed  by  the  following  series: 

oo 


n  n 


In  the  interval  in  which  stationary  state  assumption  is  justified,  |gj«l.  Hence: 

r 


1-20 


(41) 


(42) 


8 


and  according  to  Equation  17: 


i  in 


(43) 


Finally,  within  the  same  interval  {f«l9  and  Equation  40  reduces  to: 


dg 

dx 


x 

"T  “ 

-  Zee  x 


(44) 


This  new  approximate  equation  can  be  solved,  if  at  least  one  fixed  point  (boundary 
condition)  is  known.  Unfortuately,  the  initial  conditions  cannot  be  applied  for  this  case. 
Since  the  initial  value  of  g  is  g(xo)=  1,  this  function  cannot  be  considered  small  for 
x  =  Xq.  Thus,  at  the  very  beginning  of  the  polymerization,  the  approximation,  Equation 
44,  is  not  justified.  However,  it  will  be  shown  in  the  following  that  this  formal  objection 
can  be  disregarded  if  special  conditions  do  exist,  and  that  these  special  conditions  are 
also  identical  to  the  stationary  state  conditions.  Therefore,  the  following  method  will 
be  applied.  At  first  a  new  differential  equation  will  be  solved  in  which  the  function  g(x) 
in  Equation  44  will  be  replaced  by  another  function  y(x): 


dy 

dx 


*  -•£- 
T  2 

-  Zee  x  y  + 


(45) 


where  y(x0)  =  1. 

Integration  of  Equation  45  leads  to  the  following  expression: 


where 


y(x)=  -L&F{x)J  eHz)(_5__|)dz  +ffF(x0)-F(x) 
*o 

/ 


F(x)  =  2c 


x  P 
'T  "T 
e  x  dx 


(46) 

(47) 


Since  F(x)  itself  is  a  transcendental  function,  the  general  integral  appearing  in  Equation 
46  can  only  be  given  by  approximation,  or  can  be  expanded  into  a  series  which  employs 
successive  integration  by  parts: 


z  P 

f  F(z )  p  I  22P 

Je  (  —  -  i)d2  =  z  (|-  -\)e 


F(z ) 


(48) 
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Then  by  denoting 


u(z)  = 


z  P 
T_T 


e  z  z 
2c 


(49) 


such  a  transformation  repeated  n  times  leads  to  the  following  expression: 


/ 


$„(.)- R„(*> 

*  k=  I 


(50) 


where 


Sk(z)  =  u  (z)bkPk(z) 


(51) 


Pk(z)  =  ^  °kiZ 
K  i=0  K1 


-i 


(52) 


Rn(z)  =  /gF(z)d[bnun(z)Pn(z)] 


(53) 


and  bfc  and  a^i  are  constants.  According  to  Equations  50  through  53  the  function  y(x) 
given  in  Equation  47  becomes 

y(x)=4fZ  Sk(x)  +  [l_4'^  s(xo)]exp  [F(xo)_F(x)]' Rn(x,xo)  (54) 

k- I  k _ I 


where  Rn(x,  xQ)  =  Rn(x)-Rn(x0).  For  the  sake  of  convenience  Equation  54  will  be 
expressed: 


y(x)  =  t  T.  (x)-R  (x,x  ) 
k=l  K 


(55) 


where 


Ti  =  ~k  S, (x)  +  [’  "  lTSi(xo)]exp[F(xo)"F(x)] 


(56) 


T  ..=  4-Sk(x)-  ^Sk(x0)exp[F(x0)- F(x)] 


_  j_ 

k*l  ”  4 


and 


A 

Rn(x’xo}  =  /  eF(Z)bnd[un(z)  Pn(z)] 


(57) 


(58) 
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For  x»l,  Pk(x)«  0>  Hence  for  large  values  of  x: 


Sk(x)  w  bk°k,ou  (x) 


On  the  other  hand,  if  x  is  small,  xkPk(x) »  a^  k>  and 

where  a^  =  bkak,  k  *  Since  xQ  is  a  small  number, 

Sk(xo)  *  °k  [  ~ x^~  1 

Finally,  for  large  values  of  x: 


.  r  pu(xn) 1  r  i 

T,  «  --^-u(x)+  [l-  __]exp  [f(x0)-F(x)] 


and 


k*i  4 


For  small  values  of  x: 


T  ~ 

«  4x 


pu(x)  r  P«<0 


Tk^|K  Tak[Jir]  ~T0k[“7r-]  exp[F(xo)-F{x)] 


k*l 


(59) 


(60) 


(61) 


(62) 


(63) 


(64) 


(65) 


Remembering  that  y(x)  and  g(x)  are  not  identical,  the  conditions  justifying  the  stationary 
state  will  be  mathematically  applied  to  the  y(x).  To  find  such  conditions  means  to  find  an 
interval  xm|n  5.  x  £  xmax’  within  which  1y(x)l  $  1/2N.  One  has  also  to  regard  the  needs 
of  the  experimenter,  that  is,  to  have  a  stationary  state  which  (1)  appears  almost  simul¬ 
taneously  with  the  beginning  of  the  reaction  (short  period  of  induction),  and  (2)  which  is 
valid  over  a  considerable  period  of  time.  The  first  condition  enables  one  to  neglect  the 
period  of  induction  and  to  evaluate  the  experimental  curves  by  assuming  that  the 
simplifying  conditions  are  valid  practically  from  the  beginning.  The  second  condition 
ensures  that  proper  measurements  can  be  carried  out  during  a  period  of  time  in  which  a 
substantial  part  of  the  initiator  underwent  polymerization.  Such  a  stable  stationary  state 
will  appear,  if  xm^n  is  very  close  to  xQ,  and  if  xmax  is  a  large  number  which  covers  a 
substantial  period  of  the  reaction. 
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The  function  y(x)  given  in  Equation  46  will  now  be  expressed  by  the  mean  value  of  the 
integral: 


,  .  I  F(£>-F(x) 
y(x)  =  e 


[  P  ,n<T£>  -*  +  *0]  +  6 


F(x0)-F(x) 


where  xQ<£<x. 

By  substituting  x  =  xQ  into  Equation  66,  the  resulting  value  of  y(x^)  is  unity.  Therefore, 

Equation  66  satisfies  the  initial  conditions  for  y(x).  A  short  induction  period  means: 
y(xmin)  =  1/2N  also  if  xmin  is  very  close  to  xq.  Such  a  condition  can  be  achieved  only 

if  the  expression  appearing  in  the  exponent  of  Equation  66,  F(x)-F(xq),  which  is  zero 
for  x  =  Xq  will  increase  rapidly  as  x  starts  to  exceed  xQ.  Therefore,  it  is  important  to 
examine  the  function  F(x)-F(xq)  closely,  to  derive  conditions  assuming  a  short  period 
of  induction. 

MATHEMATICAL  ANALYSIS 

OF  EQUATION  46  AND  THE  EXPRESSION  F(x)  -  F(xQ),  USING  EQUATION  47, 
ACCORDING  TO  THE  STATIONARY  STATE  CONDITIONS 


Let  z  =  2t,  then  the  integral 


f~  T  “7  . 
fez  dz 


transforms  into 


,  -4  ,  -t  -4 

Kt)  =  2  z  J  e  t  dt 


Successively  integrating  I(t)  by  parts  gives  rise  to  the  following  series  expansion: 


I(t)  = 


1  --*•  p 


By  substituting  t  =  ^  z  Equation  69  becomes 


I(2)  =  T^r2"p2/22iFi<'*2-T*T>* 


where 


c  .  .  .  ,  .  ox  .  o(a-H)x2  .  a(q+l)(Q+2)x3  , 

|F|  (o»b»x)  -  bl!  b(b+l)2!  b(b+l)(b+2)3! 
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|Fj(a,  b,  x)  is  known  as  the  confluent  hypergeometric  function  (Reference  7).  Since  the 
confluent  hypergeometric  function  is  absolutely  convergent  for  all  values  of  a,  b,  and  x, 
real  or  complex,  excluding  b  =  0,  -1,  -2, . . . ,  the  integral  I(z)  always  exists  except  for 
p  =  2  +  2n(n  =  l,2,...).  In  general,  p  varies  only  within  the  interval  o i  p  S  1.  Hence, 
the  integral  I(z)  is  always  defined. 

Equations  47  and  70  clearly  lead  to: 


F(x)-F(x0) 


dz 


i- 

x 


(1,2- 


P  x  , 

2  ’  2,_Xo 


(72) 


For  very  small  values  of  x,  the  confluent  hypergeometric  function  approaches  unity,  and: 

F(*>-F<x0l*  2TTJ-[«'_T-CT]  1731 

if  o  <  xQ<  x  «  1.  Under  such  conditions,  namely,  when  x  and  xQ  are  very  small  and  x  is 
close  to  xQ,  then  F(x)  -  F(xo)  can  only  be  very  large  if  c  is  sufficiently  large,  that  is: 

F(x) -  F(x0)  >  M 

where  M  is  a  large  number,  if 


4c  > 


M(2-p) 


x  -  x0 


(74) 


Therefore,  one  may  expect  a  short  induction  period  only  if  c  is  a  very  large  number. 
Hence,  the  following  mathematical  analysis  will  be  carried  out  under  the  assumption 
that  c  is  very  large.  Let  us  now  evaluate  the  function  F(x)-F(xq)  for  large  values  of  x. 

Since  the  integral 


x  __z  _  P 

2c  f  e  2  z  ~^dz  exists  for  all  values  of  x,  then: 
xo 

X  X 


00  0 

/  *  /  *  !  *  I  *  h 

X„  00  0  X. 


(75) 
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After  repeated  integration  by  parts,  1^  yields  the  following  series  expansion: 


CI 

CD 


K  Z  P 


I,  =  2c  /  e  2  z  2  dz 


P  X 


=  -4c"V2[i-4*  PiH£Hll£l+...] 


(76) 


Therefore,  if  x  »1, 


_JL  _  x 


lf  «  -  4cx  2  e  2 


(77) 


The  value  of  Ib  is  readily  available  from  tables  of  integrals  and  is  given  as 

oo  __z  _p  s  -3L 

I2  =  2c  fe  2  z  2  dz  =  2  2cr(i--|), 


(76) 


00  .  | 

where  r  (x)  =  J  e~  T  t*  dt  is  the  well  known  gamma  function.  Since 


x0«  1,  Is  may  be  approximated  according  to  Equation  73.  That  is: 


o  __z 

1^  =  2cjezzzdz 
o 


Combining  Equations  77,  78,  and  79  one  obtains: 


4c  1  —  2 
2-p  xo 


F(x)-F(x0)  =  2.  2  cr(l--|-)-4cx  2  e  2  -  xd  2 


2-£- 

2  2  crt  l- J-  ) 


(79) 


(80) 


Since  p  «  1  and  r  (  -s)  =  JiT,  Equation  80  may  be  expressed: 


F(x)-F(x0)  =  2*  c 77,  if  X  »l.  (81) 

From  the  above,  it  can  readily  be  seen  that  the  function  F(x)-F(xQ)  approaches  its 
asymptotical  value  2a/5F  very  rapidly.  This  property  of  F(x)  is  illustrated  in  Figure  1  for 
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the  case  p  =  1  and  xQ  =  0.  By  imposing  these  conditions  and  substituting  z  = 
transforms  into  the  normal  distribution  function.  That  is: 

x  _je  _£  vT  _ 

F(x)  =  2c  f  e  2  z  2  dz  =  4c  J  e  2  dt 

O  0 

vr 

=  4c  */2~tF  J  <^>  ( t )  dt  =  2cY(  x ) 
o 


where 


+"'*  yfe  e’T 


Let  us  now  consider  the  expansion  of  y(x)  as  expressed  by  Equations  56, 


y(x )  =  T,  +  T2+  t3  -  r3 


where 


T,  =  ^uUX-f-n*  [,--L U(x0)(f 


T,  =  4u2(»W-2p- *11-  4-u*(xa)('gj~f  +i)gF(x0)~F(x) 


T.  = 


8  x2 
*0 


XuVK^-y+tP  +  -£*g_JL.„ 

8  r  yz 


*o 


R,  =  »FI;ldL*lz>  P3-4^  .  Pl|£  .  i..„l 

X0  Z  2 


t3,  F(x) 

(82) 

(83) 


57  and  58: 

(84) 

(85) 

(86) 

(87) 

(88) 
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X 


13 


X 


Figure  1.  Behavior  of  the  Function  Y(x)  as  Defined  by  Equation  83 
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Let  x  and  x  .  be  defined  by  the  following  equations: 
max  mm  J  &  ^ 


max 


and 


A 


X 


» 


max 


(89) 


and 


in  which 


4«mi„  mm  2N 


x„  <  £  ,  ,  and  5  .  Is  close  to  x„ 
o  min7  mm  o 


(90) 


Since  xmax  and  xmin  are  dependent  upon  c,  the  conditions  xmax»l  and  xmin  close  to  xq 
are  easily  satisfied  by  adjusting  c. 

Since  xmjn  is  very  close  to  xQ,  Equation  90  gives  rise  to  the  following  inequality: 


Hence: 


4x, 


u(x0) «  I 


T.  <  ^U(x)(*f--I) 


+  e  p(*o> -  F(x) 


(91) 


(92) 


F  (  x  )  —  F  (x  ) 

It  can  readily  be  seen  that  for  sufficiently  large  c  the  term  e  o  decreases 

rapidly  and  may  be  neglected  if  x  »  x  •  Then: 


T.  i 


l 

2  N 


if  "J  <  x  <  x 

min  =  x  =  max 


(93) 


where  and  x  x  are  defined  by  Equations  89  and  90  respectively.  Examining  Ta 
and  T3inthe  interval  xmin  g  x  g  xmax  one  observes  that: 


and 


T  < 
2  - 


(94) 


T  <  (  — )3 
3  =  '  2N 


(95) 
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The  integral  Rg-  appearing  in  Equation  88  may  be  expressed  by  the  mean  value  theorem  as 
follows: 


R 


3 


I  r  3,  sf  P3"4p2  +  4p  p2-2p  p  ^ 
=  TLU(XH - 7 -  +  ^2  X  j 


3,  ,/p3-4p2  +  4p  ^  pz-2p  p  ,M  _F(^)  —  F{ x) 

-u  (x0)(  ^  _ - 

*0  0 


(96) 


where  xQ  s  jsx. 
Consequently: 


1  R3>  *  (  2*f  )3»  if  \nin  i  i  * 


mm 


max 


(97) 


Since  N  is  a  very  large  number,  the  terms  Ta,  Ts  and  R3  may  be  neglected  and  the 
following  relationships  arise: 


I  y(x)l  «  lij  (x)  I  < 


.-A  ^  .  a 

if  X  ._<  X  <  xmnw 
mm =  -  max 


(98) 


Recalling  that  the  validity  of  the  stationary  state  assumption  can  be  expressed  by  the 
inequality  jy(x)(  s  where  x  .  <:  x  <;x  ,  it  is  obvious  from  Equation  98  that 

the  stationary  state  exists  in  the  interval  xmin<  x  <  xmax .  Hence,  Equations  89  and  90 

define  the  lower  and  upper  bounds  of  the  interval  in  which  the  stationary  state  exists. 

Assume  the  exponential  term  in  Equation  90  is  the  leading  one,  that  is: 


4x 


u<*min  >  «  * 


F(x0)  -  F  (x) 


min 


I 

2  N 


(99) 


This  additional  condition  on  Equation  90  forces  x  .  to  approach  x  ,  thus  decreasing 

min  o 

the  induction  period.  By  using  Equation  99,  the  expression  for  xmin  is  simplified  to: 


,  P  ,  P 

-  7  •  “7 

x  .  -  x 

min  o 


(  2  -  p  )  ln2N 
4c 


(100) 


It  should  be  pointed  out  that  the  condition  in  Equation  99  is  not  a  necessary  one.  It 
means  that  the  stationary  state  may  also  exist  without  satisfying  Equation  99.  If 
Equation  99  cannot  be  satisfied,  but  the  approximation  eXmin/2  «  i  is  valid,  the 


general  Equation  90  reduces  to  a  simplified  transcendental  equation: 


8cx 

min 


-4 


exp 


[-5? 


(  x  . 
p  min 


“-TT  I 


-  X 


I 

2N 


(101) 
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It  must  be  emphasized  that  the  stationary  state  conditions  above  were  not  derived  from 
the  proper  function  g(x),  but  from  a  hypothetical  one,  y(x).  One  of  many  ways  to  show 
that  g(x)  ca  y(x)  in  the  stationary  state  interval  is  to  prove  the  existence  of  a  point  x 
such  that  xmjn=  x  %  xmax  and  g(x)  ^y(S).  It  is  obvious  that  if  such  a  point  x  exists,  then 

lg(x)|  %  in  some  interval  containing  x,  and  therefore,  one  may  neglect  the  term 

g3(x)  appearing  in  Equation  40  in  such  an  interval.  This  neglect  of  g3(x)  in  Equation  40 

leads  to  Equation  44  and  finally  to  Equation  45,  the  defining  equation  of  y(x).  Hence,  the 
derived  stationary  state  conditions  will  be  proper  for  the  function  g(x)  if,  under  the 
given  conditions,  the  following  relationships  hold: 


or 


^  xmin^  ~  xmin^ 

'y(W  ‘  9(xmm)l  <<  |fl(xmin)l 


(102) 


Now  consider  the  interval  x_.  <  x  <  x  .  .  Since  x  .  is  close  to  x  and  x  %  1,  x  is 

o  -  -  min  mm  o  o 

a  small  number,  and  the  approximation  e  x  ~  i  holds  in  this  interval.  Using  this 

approximation,  Equation  33  reads  as  follows: 


ii 

dx 


.  -p  2 
a  -  bx  r 


By  introducing  a  new  function  v,  defined  implicitly  by: 


Equation  103  transforms  into  a  new  second  order  differential  equation: 


P_  dv 
x  dx 


2  "P  „ 
-  c  x  v  =  0 


(103) 


(104) 


(105) 


where  c2  =  ab. 

This  equation  may  be  solved  readily  for  any  p,  but,  for  the  sake  of  additional  sim¬ 
plicity  only  the  case  p  =  1  will  be  considered.  Thus,  Equation  105  reads: 


.2  ,  .  2 

d  v  .  I  d  v  c  v  _ 

dx2  x  ~S7~  x 


(106) 


By  utilizing  a  new  variable  z  defined  by: 


x  = 


(107) 
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Equation  106  is  expressed  by: 


dfv  +J_dv_w  =  o  (108) 

dz2  z  dz 

which  is  the  modified  Bessel  equation  of  zero  order.  This  eguation  has  two  independent 
solutions  I^z)  and  Kq(z).  Hence, 

v(z)  =  C,  I0(z)  +  C2K0(z)  (109) 

where  Cl  and  Ca  are  arbitrary  constants,  and  IQ(z)  and  K^z)  are  modified  Bessel 
functions  of  zero  order  (Reference  8). 


Since 


and 


dIo{z) 

dz 

dK  (z) 
o 

dz 


=  I,  (z) 

=  -K,(z) 


(110) 


(  III) 


where  (z)  and  Kj  (z)  are  modified  Bessel  functions  of  the  first  order,  the  function  r 
defined  by  Equation  104  can  be  written: 

2  C,I,(z)  -C2K,(z  ) 

r  =  2b  C,l0(z)+ C2K0(z)  (ii2) 

The  stationary  state  condition  &L.  =  0  applied  to  Equation  103  yields  the  following  rs 

UA 

function: 

i 

rs=  (^-)T  if  p  =  I  (113) 


or  in  terms  of  z: 


z 

2b 


(114) 


Dividing  Equation  112  by  Equation  114  one  obtains: 

r  _  l-g  _  C,I,(z) -  C2K, (z) 
rs  "  >+8  “  C,IJz)  +  C2K0(z) 

1 

Since  the  initial  conditions  require  r(z0)=O;  where  z0  =  2cXq  2,  Cj  =  Kj(z0)  and  C3  =  ^  (z0). 
Finally: 

_r_  _  K,(z0)I,(z)-l,(z0)K,(z) 

rs  '  K,(z0)I0(z)  +  l,(z0)K^z) 
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From  Equation  91  one  obtains: 


2cx 


i 

1 

o 


»  I 


(117) 


x_ 

if  p  =  1  and  e  «  i  . 

In  the  case  of  large  arguments,  one  may  express  the  functions  L(z),  I,  (z),  K,(z)  and 
Ka(z)  by  the  following  series  (Reference  9): 


e 

f  1  +_!_  + 

( 1  )(9)  , 

(IK9H2S)  t 

JZttz 

L1  8z 

2I(8z)2 

3!  ( 8z)3 

~  ^  1 

[l- 

(3M5) 

(3X5H2I) 

•JZttz  ' 

l  8z 

2!(8z)2 

3!  (8z)3 

ire1 

f, _ L  + 

( 1  )(9) 

(l)(  9)  (25) 

■JZ-rrz 

l  8z 

2!  (8z)2 

3!(8z)» 

m  7r6'Z 

1  8z 

(3)(5)  . 

(3K5H2I)  ( 

■JZttz 

2!  (8z)2 

3 !  ( 8z ) 3 

(118) 


(119) 


(120) 


(121) 


Using  the  above  expansions  one  obtains: 
l  -  g 

l  +  g 


e  i'+-k- 

15 

-ir  i  — 3 — 

15 

ro 

(J 

N 

15 

ir  3 

L  8ZQ 

2!(8z0)2 

Jl  8z 

2U8Z)2 

j  *  rir- 

o 

2!(8zq)2 

•I1**- 

2!(8z)2+  ] 

*2‘z°[' 

15 

lfi+-L+ 

9  +•• 

.1.  Z0‘zf,  3 

15 

9  1 

1  8zq 

2!(8z0F 

Jl  8z 

2!(8z)* 

1  *  l-eTo- 

2K8Z/ 

Ie7~ 

2!(8z)2  +  J 

Neglecting  all  terms  involving  l  / z2  and  solving  for  g,  one  obtains: 


g 


+  * 


2(Z0-2) 


( 8zn  +  8z  ) 


I  + 


3 


8z„ 


j_ 
8  z 


I  zC^-z) 

Q  V 

4  z 


(123) 
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Finally,  from  the  above  it  can  readily  be  shown  that: 


y (x  .  )  -  g^  x  ,  ) 
min  min 


g(x 


min 


«  I 


(124) 


Hence,  it  is  concluded  that  at  x  ^  the  y  function  and  the  g  function  are  almost  identical 

and  that  approximating  the  g  function  by  the  y  function  in  the  interval  x  .  <  x  $  x 

mm  max 

is  justified.  Thus  the  stationary  state  conditions  derived  for  y(x)  hold  for  the  function 
g(x)  and  express  the  true  stationary  state  conditions  of  the  system  under  consideration. 

DERIVATION  OF  STATIONARY  STATE  CRITERION 

According  to  Equation  89  the  upper  bound  of  the  stationary  state  leads  to  the  following 
inequality: 


£  x 

2  ■? 

4cx  >  N  e 


(125) 


Substituting  the  kinetic  parameters  defining  x  and  c  (Equations  32  through  36  and  40) 
into  Equation  125  one  can  obtain  the  following  formula: 


8 


S  N  exp 


[ 


Ri°8lywi 

2E  +  J 
d 


(126) 


This  inequality  can  be  considered  as  a  stationary  state  criterion  for  polymerization 
under  the  current  conditions.  In  general,  the  initial  concentration,  Iq  ,  has  also  to  be 
multiplied  by  the  factor  of  efficiency  of  initiation,  f .  For  the  sake  of  simplicity,  this 
factor  was  disregarded  in  this  theoretical  investigation. 

Equation  126  will  now  be  transformed  into  a  simpler  formula.  At  first,  the  integral 
appearing  in  Equation  13  will  be  solved,  and  the  concentration  of  initiator,  T  I],  at  any 
time  will  be  expressed: 


exp  f  RJo_(kd,o-  *d?  1 

L  E  .T  J 


(127) 


By  substituting  Equation  127  into  Equation  126  one  can  obtain  another  inequality  defining 
the  stationary  state  conditions: 


(128) 


Equation  128  shows  that  the  validity  of  the  stationary  state  conditions,  within  an  error 
l/N,  is  given  by  the  value  of  a  simple  ratio  (f  fl]  kj.)/!^),  which  has  to  be  greater  than 


19 


AFML-TR-64-405 


or  equal  to  N3/64.  It  follows  that  the  stationary  state  conditions  will  always  cease  to  be 
valid  before  the  initiator  is  entirely  decomposed.  It  is  obvious  that  the  duration  of  the 
stationary  state  increases  as  the  value  of  the  initial  concentration  of  the  initiator,  I,,,  in¬ 
creases.  On  the  other  hand,  by  increasing  the  value  of  I0,  one  will  obtain  lower  values 
for  the  average  degree  of  polymerization  of  the  resulting  polymer.  Therefore,  the  higher 
the  required  degree  of  polymerization,  the  shorter  is  the  duration  of  the  stationary  state. 

From  the  general  criterion  defining  the  stationary  state  conditions  under  the  assump¬ 
tion  that  the  temperature  of  reaction  increases  at  a  constant  rate  ,  T,  Equation  126,  one 
can  derive  another  criterion  for  T  =  O,  which  defines  the  corresponding  polymerization 
at  constant  temperature.  For  this  purpose  a  function  v  will  be  introduced,  defined  as: 


v 


(129) 


By  keeping  time  constant,  and  varying  only  T,  the  quantity  (Edt)/(RT03  )  appearing  in 
Equation  129  is  now  a  constant,  and  the  newly  introduced  function,  v  ,  is  proportional 
to  T  .  Under  such  conditions  it  is  obvious  that  if  T  approaches  zero,  v  also  approaches 
zero.  By  denoting  the  exponent  appealing  in  Equation  126  as  G,  and  substituting  Equation 
129  into  Equation  126, this  exponent  can  be  expressed  by 


Since 


*±‘ 


a,. 


(130) 


(131) 


the  exponent  G  =  -5-  k  t,  if  one  maintains  a  constant  temperature  of  reaction. 
c  d,  o 

The  appropriate  criterion  for  stationary  state  derived  from  Equations  126  and  131  can  be 
given  by  the  following  inequality: 


8 


N  exp 


(132) 


where  k^  and  k^  are  constant  during  the  reaction.  Equation  132  is  identical  to  the 
corresponding  expression  derived  previously  (References  1,  2,  and  3)  for  free  radical 
polymerization  at  constant  temperature. 

The  period  of  induction  of  the  stationary  state  can  b.e  derived  from  the  inequality 
given  in  Equation  89.  The  resulting  formula  will  be  given,  in  general,  by  a  transcendental 
equation.  However,  in  order  that  this  period  be  extremely  short,  the  general 
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expression  must  be  reduced  to  an  approximate  one,  according  to  Equation  9$  and  the 
derived  relationship  is: 


In 


2N 


<&(■' 


(133) 


By  substituting  the  original  parameters  into  c,  p,  x,  and  Xq  ,  Equation  133  can  be 
transformed  into  an  approximate  formula  defining  a  short  period  of  induction: 

In  2N^(Iokt,okdto)2  f  (,34) 

An  identical  formula  could  easily  be  derived  from  a  corresponding  expression  obtained 
in  a  previous  work  (References  1,  2,  and  3)  for  polymerization  at  constant  temperature. 
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